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ABSTRACT 


A prototype management information system called Work- 
in-Process Inventory Control System (WIPICS) has been in- 
stalled at the Naval Air Rework Facility (NARF), North 
tana. WIPICS is “expected to benefit the NARE by reducing 
rework times and rework costs. By assuming that linear 
production processes may be estimated from statistical data 
accumulated by the NARF, a linear economic model is con- 
structed which predicts a required budget in man-hours, 
material cost and overhead cost categories for a specified 
production output level. Sample problems are solved and 
parametric studies are done to determine changes in the 
required budget for restrictions in man-hours, material 
cost and overhead cost. The model will be used as an aid 
in the cost-effectiveness evaluation of WEPICS by tre 


Management Systems Development Office. 
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I. INTRODUCTION 


The Naval Air Rework Facility (NARF) North Island is 
currently in the process of installing and evaluating a 
prototype of a recently developed management information 
system. The system is called the Work-in-Process Inventory 
Control System (WIPICS) and was designed and developed by 
the ROHR Corporation, Chula Vista, California. The system 
employs the latest advances in computer technology which 
include real time data file updating and information re- 
trieval with communications via a "touch-tone" telephone 
system and remote teletypewriter terminals. WIPICS contains 
data regarding individual items and their relationship to 
the next higher or lower assembly. The system provides an 
immediate audio (pre-recorded voice) response for simple 
answers and printed reports for more detailed answers. 

Each individual item or component has a unique identifica- 
tion number (register number) which is assigned and used 
for all WIPICS transactions. As an item changes location 
or status a transaction immediately updates the computer 
Buorase record and amy subsequent queries will be answered 
with current information. 

The WIPICS formulation concept began in early 1968 when 
RADM. J. Smith (then NAVAIRSYSCOMREPAC) conducted a review 
of then current management information systems. RADM Smith 


was convinced that the Uniform Automatic Data Processing 





System (UADPS) then being implemented at all Industrial 
Naval Air Stations (INAS) was rapidly being outdated by 
the tremendous advances being made in computer technology. 
Consequently NAVAIRSYSCOM activities were no longer using 
the most modern management systems available. A "task 
force" was created and the resulting study was a systems 
design proposal for the establishment of the Management 
Systems Development Office (MSDO). 

In early 1969 MSDO and the ROHR Corporation studied 
NARF North Island to determine the applicability of an 
advanced management information system currently being used 
by ROHR. The study conclusions led to the development of a 
contract request for WIPICS by MSDO in late 1969. By July 
1970 approval had been granted by SECNAV for MSDO to contract 
for the construction of a prototype WIPICS to be installed 
at NARF North Island. The contract negotiation was completed 
in January 1971. The WIPICS prototype is now installed at 
NARF North Island and is in an evaluation stage. The ex- 


pected benefits of WIPICS are three-fold: 


(1) Allow positive control of all work-in-process 
(2) Reduction of turn-around time for repairs 


(3) Improved production rates 


Other benefits would be possible if extensions are made 


ES WIPICS: 


(1) Attendance and labor force data 


(2) Association of tooling and technical data 





(3) Integrate directly with automated storage areas. 


in addition to the approval of a WIPICS prototype in 
July 1970, SECNAV placed a requirement on MSDO to evaluate 
mee bDenefitsser WEPICS in accordance with the Resources 
Conservation Program. The documentation for this cost 
justification was to be developed by MSDO with the methods 
and data subsequently audited by the Naval Audit Service. 
Using data collected and summarized by MSDO, the Operations 
Research/Administrative Sciences Department of the Naval 
Postgraduate School agreed to conduct a study to develop 
a methodology for quantitative cost-effectiveness reviews 
Lor technological changes of this type. 

A brief comparison of the UADPS Feedback system and 


WIPICS is presented in Figure 1. 


B. STATEMENT OF THE PROBLEM 

The overall project being considered by the Naval 
Postgraduate School is the rei öffent of a methodology for 
quantitative cost-effectiveness analysis of technological 
changes at a particular organization. The results will be 
specifically applicable to the evaluation of WIPICS at 
NARF North Island, but may be generalized to any organiza- 
Perot. Which produces several products and incorporates a 
Keehnological change either in a production process or the 
control of a production process. The results of this 
cost-effectiveness analysis will determine the future of 
WIPICS as a management tool for possible use at other INAS 
as well as continued use at NARF North Island. 


10 
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The general approach to the development of a methodology 
will be to estimate production functions of all major pro- 
grams at NARF North Island both before and after WIPICS 
installation. A production function involves determining 
a specific relationship Be a selected output measure 
and selected input variables. By estimating production 
Runetlions for entire programs it is possible to avoid the 
measurement of precise relationships between WIPICS and the 
NARF North Island shops.* Two different types of production 


functions are to be estimated: 


(1) The Cobb-Douglas Production function 
de) “A Constant Elasticity Substitution Production 


function. 


To formulate continuous production functions of this type 
requires a substantial amount of raw data aggregation and 
sophisticated regression analysis techniques. In order to 
validate the models using the above techniques a third 
model will be constructed. This model will be based on 
Linear Programming Techniques and will yield a linear 
approximation tovythe production functions. 

Once the before and after WIPICS production functions 
have been estimated, the results will be analyzed to deter- 
mine the effects of WIPICS. The WIPICS effects on production 
efficiency will be translated into dollar savings for compar- 


ison with WIPICS costs. 


12 





Ger SCOPE OF THIS THESIS 

Linear Programming Techniques will be used to formulate 
two Linear Economic Models for production at NARF North 
jslanad. The two models will be for the engine program and 
the aircraft program. Both models will be examined in 
detail by using sample production requirements and conduct- 
ing parametric studies to determine the sensitivity of the 
solutions obtained. No attempt will be made to compare 
these models with the continuous models (which are being 


derived separately) in this thesis. 


15 





II. DATA ANALYSIS AND REDUCTION 


A. DATA AVAILABLE 


The Naval Air Rework Facility North Island and the 


Management Systems Development Office provided actual his- 


torical data for use in constructing and testing the linear 


economic model. 
submission of reports to higher authority and for budgetary/ 


accounting purposes. 


The data is collected by the NARF for 


For each engine or aircraft that is 


reworked by NARF North Island certain statistical informa- 


tion is accumulated: 


(1) 
(2) 
(3) 
(4) 
(5) 
(6) 
(7) 
(8) 
(9) 
(10) 
(11) 
(12) 


Type of engine or aircraft 
identificavion number 

Type of work done 

Induction date 

Production date 

Production load norm (man-hours) 
Airframe change man-hours 
Direct labor Nours expended 
Direct labor cost 

Direct material cost 
Applied overhead cost 


Navy Industrial Fund (NIF) rate 


(NORM) 


(DMHR ) 
(DLB$) 
(DML$) 
(DOH$) 
(NIFR) 


The manner in which the above statistics are accumulated 


and used by NARF North Island are enumerated in Ref. 10. 
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The MSDO personnel furnished a summary listing of the above 
statistical information for each engine and aircraft re- 
worked at the NARF for the period February 1970 through 
July 1971. The basic data consisted of approximately 1500 
observations on the engine program and 360 observations on 
the aircraft program. Each "observation" consists of a 
complete recorá of the above statistical data for a com- 
pleted NARF job. This data is listed for reference purposes 
in Appendix A of Ref. 11. In order to keypunch the basic 
data information it was necessary to assign codes for dif- 
ferent types of engines/aircraft and the types of work 


done. These codes are listed in Tables I and II. 


B. SELECTION OF INPUT/OUTPUT VARIABLES 

From the basic data it was necessary to select input/ 
output variables to be used for the linear economic model. 
A single variable consisting of NORM was selected for an 
output measure. 

Input variables were selected to be direct man hours 
expended (DMHR) direct material cost (DML$), applied over- 
head cost (DOH$) and number of days in shop (NDAY). The 
number of days in shop was defined to be production date 
minus induction date. These input variables were selected 
because they were thought to represent the variable costs 
Siaproauction. | 

The expected benefits of WIPICS as stated in Chapter I 


would be expected to cause measurable decreases in the 
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Table I 


Codes Assigned to the Engine Program 


Eupine Type 
T58-GE-1 (A/F) 
T58-GE-3 (A/F) 
T58-GE-5 (A/F) 
T56-A-8P 
T58-GE-5 (C/G) 


Ј57-Р-ҢА/22 
J57-P-10 
J57-P-20A 
J57-P-22 
J57-P-420 


J79-GE-8B 
J79-GE-8B/C 
J79-GE-8B (RDTE) 


J79-GE-10 


Weck Type * 
OVHL 
OVHL/CONV 
PAR/REP 
SUP/REP 
SUP/REP/CONV 


Code 





51 
99 
55 
56 
57 


al 
T2 
T3 
74 
i 


91 
92 
JS 
95 


Code 
01 
02 
03 
04 
05 


Engine Type Code 
T6O4-GE-3 (A/F) 63 
T64-GE-6B 65 


T64-GE-6B (PAWN) 66 


T64-GE-7 (A/F) 67 
T64-GE-413 69 
T58-GE-8F 81 
T58-GE-8B 82 
T58-GE-8B/F 83 
T58-GE-8B (C/G) 84 
T58-GE-8B (HH2) 85 


T58-GE-8B/F (CONV) 86 


T58-GE-10 89 


* Abbreviations are defined 


as follows: 


OVHL 
CONV 
PAR 
SUP 
REP 
SEA 


CONUS 
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overhaul 

conversion 

planned aircraft repair 
supply 

герат 

South East Asia 


Continental United 
States 





Table II 


Codes Assigned to the Aircraft Program 


Aircraft Type Code Aircraílt Type Code 

C-2A 10 CH-3B 31 

E-2A/B 11 RH-3A 32 
SH-3A 33 

F-4J pu SH-3A/G 34 

F-4B 22 SH-3D 35 

F-4G/B 23 CH-46A 41 

F-8J 25 CH-46D 42 

F-8H 26 CH-46F 43 

RF-8G 27 UH-46A 4 4 
UH-46D 45 
CH-53A 48 

Work Type** Code CH-53D 4g 

OVHL 1% 

PAR 2% 

PAR/CONV S 

PAR/MOD 4# 

PAR/MOD/REP 25 

PAR/REP 6* 

PAR/SEA ps 

PAR/CONUS 8% 


* Second digit refers to the cycle number of the aircraft. 
(A zero means not applicable.) 


** See Table I for Work Type abbreviations. 


D 





input variables per unit of output demanded. These de- 
ereases can in turn be associated with dollar savings and 
therefore used as a measure of effectiveness for WIPICS. 

As a first step in determining the relationships between 
the available raw data elements, a computer program was 
written to compute the correlation coefficients between 
each pair of variables. The results are shown as Tables III 


and IV. The computer program is described in Appendix A. 


* 


C. CLUSTER ANALYSIS 

To formulate a linear economie model it is necessary to 
define a process. A process consists of a mathematical 
description of the amount of resources used (i.e. labor, 
material, etc.) to produce a unit amount of an output meas- 
ure. The problem is how are the raw data observations used 
to estimate a process. The most natural way would be to 
estimate a process as an average amount of resources re- 
quired to do a specific type of work on a certain engine 
or aircraft. This procedure may not be desirable because 
there may be alternate processes available to produce the 
same output. For example, to overhaul an F-8J it may be 
possible to use a "normal" amount of man-hours and have a 
low material cost. An alternate process may exist whereby 
the same work may be accomplished by spending more money 
on material and reducing the man-hours required. 

In an attempt to identify these alternate methods within 


the basic data observations a clustering algorithm was 


18 





Ceoco" T 
80850*0 
¥27SS °0 
28665 0 
Eb 26 0 
91856 0 
65819*0 
AVON 


BP4eq M?H eu[Sug doj SYUSTITFIJSOH цотзетәааод 


80850 *0 
00000°T 


00260 °0- 


ECEVS°O 


YLEOT? 0- 
CT90 T = 
9esvo*o- 


UdIN 


9225S °C 


00260 *0- 


00000°T 
88269°0 
60686 *0 
cEC66°0 
v98£6*0 
$HOQ 


S3N19N3 


L86S9°C 
Ec ESSO 
88¢S9°0 
00000°T 
69189 *0 
821659"*C 
0c979*0 
$ WO 


с82%6° 0 


vHSGONEO- 


60686*0 
69189 *0 
00000°T 
S6966°0 
67696 °0 
$910 


"ІІІ КОКС, 
19° OT 90° ЕЕ 
O€*IZOT — CO*GISZ 
9L'll62 06°5%0% 
LL°96€S  €0**91G6 
22 "8622  T9*OEZE 
T6°9LE 2т° 55 
6S LIE SZ*vES 
AJQ OLS МУЗИ 
9184S°0 6%8L4°0 
01901°0- 925%0*0- 
260660 %98Е6° 0 
821590 02919*0 
G6966°0  $%646°0 
000C0°T  z€8%6°0 
ZEB%6°0  00000°T 
YHWO WYON 


- VIVO МУУ 303 NOILVI3YWAO)I 


AVON 
YA IN 
$HOQ 
$1WQ 
$810 
YHWO 
WYON 


AVON 
YAIN 
$HOQ 
$ 1W0 
$810 
YHWO 
WYON 


L9 





00000°T 
88862°0- 
CORB "C 
Gl6LS°0 
86218°0 
L2518°0 
SS098°0 
AVON 


v?jeq MeM 3JedOoJ[V JdOJ S3UeTOTJJe0O9 UuoTAa?TeddO) 


88862*C- 
00000°T 
15126 *0- 
6598€*0 
6619€ *0- 
9€99€*0- 
¿ ¿9 ¿2° 0= 
*dIN 


¿c€6y%8° 0 
1ST2€°0- 
00000 ° I 
I7019°0 
£c696 *0 
с9516°0 
10856 0 
$HOQ 


LAVYIAIV 


$161S*0 
6598€£”*0 
1701970 
00000*I 
99119*0 
€6T29°O 
66109 °0 
$ IWA 


86218*0 
6619€°0- 
€2696°0 
99LT9°O 
00000°T 
G2766°0 
S8S96°0 
$810 


oe. Ve 
20 ZEI 
ev "81812 
EI-SOCTI 


66 *I824c 


l6*S8t* 
eS “7047 
A30 01S 


12518°0 
AA 
296L6*0 
6729 `0 
€€966°0 
oco00”T 
€cTlL6°O 
YHWOQ 


-= VIVO МУУ 3043 NOILVI3UAOI 


"AI ƏTAEL 


BG" CU 
OT°6cST 
0021655 
6S °SSTIZ 
OL°98BCS 
92° 8558 
0S °2S€8 
NV3W 


SSC9B°O 
€coce 0 
10856* 0 
66109*0 
S8S96° 0 
DOLES 
00000 °T 
WYON 


AVON 
Yd IN 
$HOQ 
$1WQ 
$910 
YHWO 
WWUON 


AVON 
Y3IN 
$HOQ 
$ 1WQ 
$870 
YHWO 
WYON 


20 





developed based on a method described in Ref. 8. A series 
of four computer programs are necessary to achieve the 
final results. These programs and the purpose of each are 
described in Appendices B through E. The general procedure 
is to select a group of observations containing several 
different types of engines (or aircraft) and some different 
types of work within each engine group. These observations 
are then analyzed mathematically for similarity of a 
selected vector of input/output measures. Similarity 
among observations is defined as those vectors that use 
input resources in approximately the same fixed proportion 
mer a unit output measure. 

Several trial selections were made with both engine 
ema aircraft data. The engine data was observed to cluster 
groups of the same engine type together with a high degree 
of consistency. For example if 20 J-79 engines were 
selected as part of a trial run of 100 observations, a 
typical result would be for 17 or 18 of the original 20 
LO be clustered in the same group. However, the trial 
selection with the aircraft data were not nearly so con- 
sistent. The cluster analysis showed the aircraft data 
had wide variations among the input variables. The clus- 
terings observed were highly irregular and containing a 
Mixture of different aircraft and different work. 

Based on the results of the engine cluster analysis a 


decision was made to form processes consisting of only 


er 





similar engines (or similar aircraft) with identical types 


of work. 


D. DATA AGGREGATION 

For the engine data the identification number consisted 
of the job order number used at the NARF. The job order 
number consists of a series of alphameric codes as specified 
in Ref. 10. Using these codes the entire engine raw data 
deck was sorted into job order number. Similar engine/ 
work types were then grouped together within each job order. 
The resulting raw data deck was then ordered by engine 
type, job type and calender quarter in widen tHe work was 
done (in special cases where a very few observations were 
available the engine/job types were not separated into 
calender quarters). Each calender quarter then represented 
a separate "process observation" on the NARF. At this 
point the data deck showed 47 different engine/work types 
and 103 processes by which the work could be done. A 
computer program was used to examine the dominance rela- 
tionships among all processes with the same engine/work 
type. This computer program is described in Appendix G. 
The dominance relationships among processes are very impor- 
tant in a linear program. For example, if process 1 and 
poocess 2 produce the same output but process 2 uses more 
of every required resource than process 1, then process 1 
dominates process 2. In a linear program that minimizes 


costs, process 1 will always be selected as preferred to 
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process 2. Therefore, if the linear programming model is 
to have alternate processes that reflect choices to be 
made then the dominance must be eliminated between all 
processes with the same output. 

The computer program in Appendix G was run with the 
103 processes as previously determined based on job order 
numbers. The dominance relationship between each pair of 
processes that produced the same output was examined. If 
no dominance existed then that group of alternate processes 
was left alone. If dominance existed among any of the 
processes in the group then a subjective decision was made 
to combine two of the processes. For example, if there 
were five processes and dominance existed between three 
pairs of them, then two processes would be selected to be 
combined. The decision on which pair to combine was 
entirely subjective. The group would then have four alter- 
nate processes remaining. The dominance computer program 
would then be run again and the en alternate proc- 
esses would be examined for dominance. This procedure 
was repeated until all dominance was eliminated within 
process groups that had the same output. Since the outputs 
were intended to be measured separately (as opposed to an 
aggregate output measure) it was not necessary to examine 
the dominance between processes of different outputs. 

In combining individual raw data observations into 
estimated processes an aggregation computer program was 


used. This program is listed in Appendix F. The program 
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estimates the process by calculating an average value for 
each variable on the raw data cards. The resulting engine 
raw data was aggregated from 103 processes to 81 processes 
in order to eliminate all "within" group dominance. 

The aircraft data consisted of 365 raw data observations. 
Fach type of aircraft and type of work considered as a 
separate output meant that 70 outputs were required. This 
small amount of data then had to be used to estimate not 
less than 70 distinct processes. A decision was made to 
estimate only a single process for each output because of 
mae limited data available. In this case there were no 
alternate processes so it was not necessary to check for 
dominance. The linear program in this case would not have 
a choice among processes and therefore the solution becomes 
Bela, Further discussion of the aircraft data is con- 
tained in Chapters IV and V. 

The final listings of the estimated processes for the 
aircraft and engine programs are оаа as Appendixes H 


and I respectively. 
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III. LINEAR ECONOMIC MODELS 


A. INTRODUCTION 

Linear programming is a mathematical technique for 
solving constrained optimization problems of a special type. 
A linear economie model is a specific type of linear program 
consisting of the maximization (or minimization) of a line- 
wr unction of n variables subject to m linear Inequality 
constraints. Most unconstrained optimization problems may 
be solved by the determination of critical points and their 
nature (maximum, minimum or saddle point) through solution 
of first and second derivative equations. In the case of 
EXE"Inear funetionothe first derivative is a non-zero constant 
and all higher order derivatives are zero. The addition of 
a constraint set to the linear function further complicates 
Mathematics. However, several mathematical tools have 
been developed to solve optimization problems where the 
objective function and the constraints are linear functions. 
The first explicit use of linear programming techniques 
weeurred in connection with planning activities for the 
U.S. Air Force in the late 1940's. Linear programs (and 
their close mathematical relative the input-output models) 
have become increasingly important to microeconomic theory 
im recent years. 

In the linear economic model being constructed it is 


necessary to define a process. A linear production process 
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is an activity by which one or more outputs are produced 

im faxed proportions» by the application of one or more 
раев in fixed proportions. The production process as 
defined is homogeneous of degree one which implies constant 
returns to scale. For example, if all inputs to a process 
are doubled then the output will also be doubled. A linear 
meoauction function ís formed from a collection of linear 
production processes that may be used simultaneously. 

The optimal solution to a linear economic model as 
described above consists of finding the combination of m 
processes from n available processes such that a linear 
objective function is maximized (or minimized). Detailed 
information concerning linear programming and microeconomic 


Dheory can be found in Refs. 1 through 5 and in Ref. 9. 


B. ASSUMPTIONS 
The following assumptions are made in the formulation 


of a linear economic model: 


(1) The estimated processes (reworking of engines and 
aircraft) are linear functions and therefore these processes 
exhibit constant returns to scale. 

(2) The above linear processes may be estimated by the 
aggregation of a finite set of observations over some time 
period. 

(3) The NARF is not a profit maximizing organization, 
therefore the management objective of the NARF will be 
assumed to be minimization of costs subject to completion 


of all work demanded by the operational forces of the Navy. 
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(4) Prices used in the model are constant and may be 


estimated from the production data furnished by the NARF. 


C. FORMULATION OF THE LINEAR ECONOMIC MODEL 
To present the model it is necessary to have some precise 
mathematical notation available. Definitions will be as 


follows: 


m number of output constraints 


n = number of production processes 
(n 2 m) R1 
R = column vector of available resources = |R2 
R3 
where 


Rl = resource for direct man-hours expended per man-hour 
of NORM 

R2 = resource for direct material cost in dollars per 
man-hour of NORM 

R3 = resource for direct overhead cost in dollars per 
man-hour of NORM 

P = column vector of prices associated 


with Rl, R2 and R3 respectively = 11.0 


where p = price of labor in dollars per man-hour 


Z = column vector of activity levels = 
des len 
To teenucdYopy matrix of “observed 
processes at NARF = | 


[(m+3) x n matrix] 


e 





where C; 


ПОШ ОШО ШК 
units of output (NORM) 


01. 150 N 
(mxn matrix) O ODO 
0 OL 


resource matrix = [715] 
I= l, 2,3 
ЕЗЕП 
with rij = amount of quen resource used per unit 


outpHutoob m activity. 


column vector of production 
output desired in man-hours 
of NORM. = [Y:] 


jl = lm 


diagonal matrix of number 


T. 
of days in shop per unit rt O 


of output = © E 
Un 


column vector of penalty 
costs = [ey] 
a уп 


—Mpenaloy cose in dollars incurred for “each day an 


shop associated with the ¡Eh process. 
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Now that a basic mathematical vocabulary exists for the 


variables, processes and costs, the linear economic model 


will be presented and discussed: 


ebjective function minimize PIT, m Cm, 
constraints -T, _y 
subject to: |------- ° 7 qL nei acm 
T, R 


Ehe objective function consists of two terms: 


p oT, 7 


2 2 aotusludost incdollassaofonesouPcesauscdaebos 
activity level vector Z 
с.т. = penalty cost in dollars for activity level 


vector Z. 


The constraint set may also be separated into two parts for 


discussion: 


=T, :2 < -Y is a set consisting of the first m rows ot 


the constraints which force the activity 


Tevel vector 4 Lor choose a"setrotiprocesses 


Which satisfy the production vector y. 


Е 
сз 
IA 
"y 


consists of the last 3 rows of the con- 


straint set and ensures that actual resources 


used do not exceed available resources. 


D. SOLUTIONS OF THE MODEL 


It is clear that the size of the resource vector R is 


highly dependent on the production vector Y. 


In fact these 


two vectors must be chosen carefully to avoid an infeasible 
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linear program. The easiest approach to avoiding an infea- 
epee linear programms tommake the resource Vector R very 
large, therefore the optimal solution to the model will not 
me constrained by an active resource constraint. In economic 
terms this means that the "customer" is willing to pay the 
NARF as much as necessary to accomplish the work specified 
by production vector Y. If the NARF is assumed to be effi- 
ciently operated, then the costs incurred will be minimized. 

In summary, the objective function of the model will be 
a minimum when the model is not constrained by the resource 
Mector R. In other words, the costs are minimum when the 
fmeseurce constraints are inactive constraints. The value 
of the objective function represents the total cost incurred 
to accomplish production vector Y based on past performance 
Gata of the NARF. 

It should be noted that the linear economic model pre- 
sented in this chapter is not a production management tool 
in the sense that the model chooses "processes" by which the 
NARF should rework aircraft or engines. The model only 
puvides a budget cost plus a penalty cost for a specified 
amount of work to be done. Some information is available 
from the model concerning tradeoff values among the three 
resources. The model can also determine the minimum amount 
of any single resource required to accomplish the specified 
work. Some example solutions are presented and discussed 


in the next chapter. 
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IV. SAMPLE PROBLEMS 


A. INTRODUCTION 
Two separate models were constructed using forecasted 
workload requirements during third quarter, fiscal year 
1972. The forecasted requirements for the engine and air- 
craft programs were obtained from Naval Air Rework Facility 
Production and Planning Notices dated 17 June 1971 and 26 
July 1971 respectively. Тһе forecasted requirements were 
converted into two separate production vectors (Y's). 
Penalty costs were computed for each type of engine 
and aircraft reworked at NARF North Island using information 
supplied by MSDO. Aircraft penalty costs were computed by 
ppopating the total cost of the aircraft over its estimated 
service life. The resulting penalty cost was in units cf 
dollars per day in shop. A similar computation was desired 
for engine penalty costs. However, an estimated service 
Tife for each engine was not available. When engines 
are reworked by the NARF, engine components are replaced 
as necessary resulting in a "new" engine. For this reason 
the service life of all engines was assumed to be ten years. 
Table V shows the computed penalty costs. For convenience 
Tables V through XI are located at the end of this chapter. 
A constant was also needed for converting direct man 
Pours into a dollar cost. As previously defined in Chapter 


ПСА Che value of p in the vector of prices, b. 
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A computer program was developed to compute the average 
labor cost. The computer program is shown in Appendix J. 


The results were as follows: 


._P (in dollars per man-hour) 
engines 6.04 


aircraft 6:17 


B. COMPUTATIONAL PROCEDURES 

To solve the engine and aircraft linear programs Mathe- 
matical Programming System/360 (MPS/360) was used. This 
IBM applications program is composed of a very complex 
set of computational procedures. A subset of these proce- 
dures are used to solve linear programs. However, the user 
meee Construct a control program that orders a specific 
sequence of procedures to be executed by the computer. A 
detailed knowledge of Refs. 6 and 7 is necessary to use 
MPS/360. 

In addition to preparing a control program for MPS/360, 
the linear program input data must be in a special format 
described in Chapter 4 of Ref. 7. To assist in preparing 
data in this format a computer program was prepared. This 
program takes specified input data and produces a deck of 
punched cards for MPS/360 in the required format. This 


program is documented in Appendix K. 
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C. ENGINE MODEL 
1. Input Data 

The sample data used for the engine model is listed 
u Rable VI. Each row of the T, submatrix of matrix T is 
associated with a specific type of engíne and type of work 
done on that engine. The production vector Y is associated 
with each row of the T4 matrix because Y specifies the re- 
quired number of man-hours of NORM for the particular type 
of work. The linear program must then choose the process 
(columns of matrix T) by which the cost of the work will 
be a minimum. If only one process is available and produc- 
tion is required then the linear program must select that 
process. 

тре toval amount of resources available are shem 
as rows R1, R2 and R3. These values form the resource 
Мес гор К. 

The model was solved for three different, but 
су related, objective functions. The values listed In 


columns OBJ1, OBJ2 and OBJ3 of Table VI are calculated as 


follows: 
OBJ1 = Р.Т, 5 7 = total cost per unit output of 
NORM (total cost = resource cost + penalty cost) 
OBJ2 = PTT, = resource cost per unit output of NORM 
OBJ3 = cT, = penalty cost per unit output of NORM 


In Table VI OBJ1 


OBJ2 + OBJ3 except in cases where numbers 


Meme rounded Oli by the computer. 
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Recall from Chapter TIT that the T; matrix consisted 


of a 1 in only one row of each column of T, and all the 


т 
other elements of that column were zeros. The row in which 
the 1 belongs is given in Table VI in the T1 ROW column. 

The columns in Table VI labeled Rl, R2 and R3 are the three 


elements of the T, matrix described in Chapter III. These 


2 
three values correspond to the amount of each resource 


Mequired per unit output of NORM for that particular process. 


The column labeled T3 in Table VI is the number of days in 


shop per unit output of NORM for that particular process. 


The values for vector C are found in Table V. However, 
the penalty costs per unit output of NORM may be found 
directly in the OBJ3 column of Table VI. Recall that 
OBJ = cl. 
3 
BRO. Solutions 

The solutions to the three linear programs (each 
with a different objective function) using the engine model 
data are presented in Tables VII, VIII and IX. The optimal 
solution for each problem is shown in the first column of 
each Table. The actual choice of processes selected by 
MPS/360 is not shown, because as stated in Chapter III, 
these "processes" do not represent actual physical procedures 
which are used by the NARF. The processes do give an estimate 
of the total production cost for a specified production 
vector. The anticipated change in production costs will be 


Ghe foundation for the cost justification of WIPICS. 
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It is necessary to explain the relationship between 
Ehe Optimal solutions presented in Tables VII and VIII. 

The solutions are nearly identical. "Mis result isnot 
unexpected as can be seen in Table VI. The magnitude of 
OBJ2 (resource cost) is much larger than OBJ 3 (penalty 
cost) in every case. By attempting to minimize the total 
cost of production (resource cost + penalty cost) the 
linear program arrived at a solution that also happened to 
omen ze resource cost. This is not an uncommon result 
when one component of an objective function is much larger 
than another component. One might adde observe that the 
value of the objective function in Table VII is approximately 
90% resource cost and 10% penalty cost. This is another 
reason for the two solutions to be similar. 

In Table IX the solution to minimization of penalty 
cost turns out to be very costly in terms of required re- 
source costs. An increase of nearly $300,000 in resource 
cost for a penalty cost reduction of only $20,000 is the 
difference between Tables VII and IX. The obvious conclu- 
sion is that unless penalty cost is higher than presently 
psuNculated, the resource cost will nearly always determine 
the optimal solution. A discussion of alternate methods of 
determining penalty cost will be discussed in Chapter V. 

3. Sensitivity Analysis 

Using procedures described in the MPS/360 references 

a parametric study was done on the resources DMHR, DML$ 


and DOH$. Since the optimal solution was obtained with 


35 





Uneonstrained resources available it is not necessary to 
increase any resource. The optimal solution shown in Tables 
VII, VIII and IX will remain the same for all values of 
resources greater than those shown in Table VI. It is pos- 
sible to obtain different solutions to each model by de- 
creasing the resources one at a time until the linear program 
Becomes infeasible. Infeasible means that there is not 
enough of some resource in R to accomplish the work specified 
by the production vector Y. In other words, a confliet 
exists between some of the constraints in the model and no 
solution exists for that linear program. The results of 
decreasing each resource (one at a time) are shown in the 
parametric studies section of each solution table. 

If the available quantity of a resource is reduced 
in value until the entire quantity is consumed in production, 
then any further reduction below that value will result in 
an increased cost in the objective function. In linear 
programming terms the shadow price of the dual variable 
associated with that resource constraint is frequently 
positive when the constraint is binding. The increased 
cost of production as described above can be observed in 
maples VII through IX. 

An interesting result appears in Tables VII and 
VIII. When the DMHR resource was reduced, the DOH$ resource 
was also reduced in the new solution (but the DML$ resource 
use increased). The explanation for this effect is found in 


Table III. The correlation coefficient between DMHR and 
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DOH$ is +.99032. This high correlation suggests that over- 
head costs at the NARF are applied to each job based on the 
DMHR for the job. This is the way overhead costs are 
applied at NARF, North Island. Therefore, a decrease in 
available DMHR could also cause a reduction in DOH$ used. 
In summary, the parametric studies on the resource 
Vector indicate the minimum amount of available resources 
mequired for a feasible linear program. However, if all 
three resources are reduced to their respective minimum 
values at the same time an infeasible linear program will 
result. The reason for this may be seen in Table VII. In 
reducing DMHR the requirement for DML$ is increased. The 
parametric study results as shown were used to constrain 
only a single resource and the two remaining resources were 
considered available in unlimited amounts. Other types of 
parametric studies are available if desired. For example, 
mm hree resources may be changed in fixed proportions. 
The MPS/360 Users Manual (Ref. 7) discusses available post- 


optimal procedures. 


D. AIRCRAFT MODEL 
l. Input Data 
The sample data for the aircraft model is very 
similar to the engine data described previously. Each row 
of T4 is associated with a particular type of aircraft and 
type of work. The procuction vector Y and the resource 


vector R are the same as described in the engine model. The 


data is listed in Table X. 


37 





The aircraft moder consisted of a single objective 
function which was calculated as follows: 


OBJ = m Т, + ВЕ Т; = total cost per unit output of NORM 


(total cost = resource cost + penalty cost) 

The remaining columns šn Table X are identical in 

meaning to the engine model descriptions for Table VI. 
СУ Solution 

The solution of the aircraft model is presented in 
Table XI. Due to the manner in which the data was aggregated 
me aircraft model consisted of only one process for each 
ne of Job in the production vector Y. The T; matrix 
consisted of 70 rows and 70 columns. Linear programming is 
not necessary to solve a problem of this type. There are no 
choices to make among processes, therefore, the activity 
level vector Z must be identical to the production vector Y 


(because T, is an identity matrix). The solution in Table 


1 
XI does give an estimated total production cost associated 
Eh the production veetor Y which can be used for budget 
estimates. A discussion concerning the aggregation problems 
for the aircraft data is presented in Chapter V. No sensi- 
tivity analysis was conducted on the aircraft model because 
the amounts of resources required at the optimal (uncon- 
strained resources) solution also represent the minimum 
resources required for a feasible solution to the linear 


program. This result occurs because there is only one process 


Bor cach job type in the production vector Y. 
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Table V 


Penalty Costs 


Pircraft Program 


Type Unit Service Life Penalty Costs 
Aircraft Cost (months) ($/day in shop) 
C-2A $10,742,000 105 280.29 
E-2B 10,742,000 105 280.29 
F-lB 2,756,000 80 94.38 
RF-4B 2,699,000 96 77.02 
F-4J 2,492,000 80 , 85.34 
RF-8G 1,300,000 15 47.48 
F-8H 1,352,000 54 68.59 
F-8J 1,302,000 54 66.05 
RH-3A 1,508,000 80 51.64 
SH-3A 19279006 80 43.76 
SH-3D 1,064,000 80 36.43 
CH-46A 1,083,000 80 37.08 
UH-46A 834,000 80 28,56 
CH-46D 1,083,000 80 37.08 
UH-46D 834,000 80 28.56 
CH-46F 1,083,000 80 37.08 
CH-53A 1,453,000 96 41.46 
CH-53D 1,453,000 96 41.46 
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Table V - Penalty Costs (continued) 


Enpine Program * 


Type Unit Penalty Costs 

Engine Cost ($/day in shop) 
J79-GE-8B $225,000 | 61.64 
J79-GE-8C 225,000 61.64 
J79-GE-10 180,000 49.32 
T58-GE-8B 65,000 175817 
T58-GE-8F 65,000 JO 
T58-GE-10 65,000 Ligon 
T64-GE-6B 125 ,000 34.25 
T6li-GE-113 167,000 45.25 
J57 (all models) 210,000 51.53 
T56 (all models) 98,600 27:01 


* a 10 year service life was assumed for all engines. 
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PRODUCTION VECTOR ( Y ) 


TYPE/WORK CIDE 


ROW NO. 
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Table VI. Engine Model Input Data 
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RESOURCE VECTOR ( R ) 
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Engine Model Input Data (Continued) 


Table VI. 
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Engine Model Input Data (Continued) 


Table VI. 
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PRODUCTION VECTOR ( Y ) 


TYPE/WORK CODE 
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Aircraft Model Input Data 


Table X. 
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RESOURCE VECTOR ( R ) 
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Aircraft Model Input Data (Continued) 


Table X. 
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Aircraft Model Input Data (Continued) 


Table X. 
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Aircraft Model Input Data (Continued) 


Table X. 





Table XI 
Aircraft Model Solution 


for Total Cost Objective Function 


Optimal Solution with Unconstrained Resources 


Objective Function Value 12,400,594 
Resource Cost 11.69] 157 
Penalty Cost 709,443 
DMHR Used cS 
DML$ Used 1,921,528 
DOH$ Used 4,983,956 
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V. CONCLUSIONS AND AREAS FOR FUTURE STUDY 


A. CONCLUSIONS 

The sample problem solutions in Chapter IV answered 
some basic questions concerning the operation of a NARF. 

In the beginning of this study, it was not known if trade- 
offs existed between resources used by the NARF. The solu- 
tions in Table VII showed that the ranges of the resources 
used were rather limited for feasible solutions. For exam- 
ple, in Table VII the reduction in use of DMHR to a minimum 
resulted in a 2.3% decrease from the optimal unconstrained 
solution. Similarly minimizing DML$ and DOH$ produced 
resource reductions of 4.2% and 11.2% respectively. Asso- 
ciated with minimizing DMHR, DML$ and DOH$ (one at a time) 
the total cost of production increased approximately 2 to 
3% in each case. This demonstrated that the available 
tradeoffs among resources are small in the linear economic 
model. This result was not unexpected because of the simi- 
larity among processes estimated for each type of work. 

The total cost objective function is composed of two 
components, resource cost and penalty cost. In the sample 
problems considered the total cost was approximately 90% 
resource cost and 10% penalty cost. This fact caused the 
total cost optimal solution to be the same as the resource 
cost optimal solution. Alternate methods of determining 


penalty costs are discussed in Section B of this chapter. 


De 





The sample problems solved in Chapter IV predicted 
budget requirements for a specific calender quarter, based 
on forecasted production requirements. To draw a a 
ful conclusion about the validity of the model these results 
must be compared with the actual production work in that 
quarter and the actual quantities of resources consumed. 
This type of comparison should be made over a series of 
production periods of varying length to determine the 
accuracy of the model. This validation was not completed 
uum ime for inclusion in this thesis. Since this modelis 
only a part of a continuing project at the Naval Post- 
graduate School, additional results are expected in the near 


future. 


B. AREAS FOR FUTURE STUDY 

1. Expanding Number of Observations 

As mentioned in Chapter II and IV, special problems 

were encountered with the aircraft model sample data because 
of the small number of raw data observations. Based on 
discussions with NARF and MSDO personnel a decision was 
made to consider work done in different "cycles" (a cycle 
is the number of times the aircraft has been through the 
NARF for the same type of work) as different types of work. 
The resulting aircraft model had 70 separate processes 
with only 365 observations to estimate the processes. In 
comparison the engine model had 81 separate processes and 
1532 observations. Forty-one out of 70 processes in the 


aircraft model had three or fewer observations. The engine 
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program had only 11 out of 81 processes with three or fewer 
observations. The magnitude of this problem could be re- 
duced by increasing the number of observations in each model. 
However, an upper limit should be placed on the age of an 
observation to ensure the processes that are estimated 
reflect the current technology at the NARF. 
2 Penalty Costs 

As SE in Chapter IV the penalty costs had very 
little effect in determining the solution to the engine 
model. The reason was that the resource cost was at least 
five times as large as the penalty cost when the objective 
funetion was formulated. (See Table VI for examples). 
It could easily be argued that the present method used to 
determine the penalty costs does not approximate the 
opportunity cost (for the Navy) of keeping an engine (or 
aircraft in the shop for an extra day. 

One alternate method by which the penalty cost 


might be computed is: 


procurement cost 
+ 


total maintenance cost over lifetime 


penalty cost = 
lifetime 


Certainly there is not a single procedure to determine 
penalty cost that cannot be criticized in some area.  How- 
ever, this problem must be addressed and resolved if an 


accurate cost-effectiveness analysis of WIPICS is to result. 
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Зе Production Vector 

In the linear economic model formulated in Chapter 
III the production vector (Y) is assumed as a given input 
to the model. In reality, the production demanded of the 
NARF each calender quarter is the result of a complex 
negotiation procedure between NAVAIRSYSCOM representatives 
and NARF representatives. Each component (Yi) of the 
production vector may be considered a random variable 
Кобе Value depends on the negotiation results. It is 
possible to construct a mathematical model of the negotia- 
tion procedures that would generate (predict) production 
ED tcpPs for use in the linear economic model. Further 
consideration should be given to a model of this type in 
lieu of using fixed production vectors based on historical 
data. 

H. Management Constraints 

One area in which extensions of the linear economic 
model would be beneficial is the incorporation of con- 
straints for management practices. For example, there is 
a fixed manpower ceiling at the NARF for the permanent 
(employed at the NARF for more than one year) labor force. 
BNEDOSse a given production vector (Y) utilizes the entire 
DMHR resource making the model infeasible. A temporary 
labor force could be hired from the surrounding community 
and be used to make up the deficit DMHR. The model becomes 


feasible. 
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Another possible extension exists in the pricing of 
labor. The management policy at the NARF could dictate 
a maximum allowable percentage of overtime per quarter. 
The overtime man-hours could be modeled as an additional 
Ponstraint and priced according to current union contract 
agreements. 

Further study is required to determine the effects 
of management decisions (or management constraints such as 
union contracts) on the linear economic model. Modifica- 


tions and/or additional constraints may be necessary. 
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APPENDIX A 


Correlation Program for Raw Data 


pose: The purpose of this program is to calculate 
and print a table of correlation coefficients between 
NORM, DMHR, DLB$, DML$, DOH$, NIFR and NDAY. A table 
of means and standard deviations is printed for each of 


the above variables. 


B. Inputs: The first card is a control card which is then 
followed by the raw data cards. All fields are right 
adjusted unless otherwise specified. The abbreviation 
"CC" refers to card columns and will be used in all 


Appendices. The format of the cards follows: 


Control Card: CC 1-4 number of raw data cards 
CC 5-8 0001 for aircraft 


0002 for engines 


Raw data card: CC 1-2 engine/aircraft code 
CC 3 blank 
CC 4-10 Identification number 


CG Ad blank 

CC 12-13 work code 

cc 14 blank 

CC 15-18 induction date 
cc 19 blank 


CC 20-23 production date 


zu 





CC 21 blank 

CC 25-29 NORM 

cC 30 blank 

CC 31-35 alrframe change man-hours (blank 

for engine cards) 

CC 36 blank 

CC 37-41 DMHR 

СС 42 blank 

cc 43-48 DLB$ 

cc 49 blank 

CC 50-55 DML$ 

CC 56 blank 

CC 57-62 DOH$ 

Cer os blank 

CC 64-68 NIFR 


The NIFR is expressed as an integer but the last two 
digits are cents. For example a-‘NIFR of 1482 means $14.82. 


All other dollar values are in whole dollars. 


Outputs: A table of correlation coefficients and a table 
of means and standard deviations for each variable are 


printed. 


special Notes: 

(1) Subroutines CORRE, LOC and MSTR are from the IBM 
Scientific Subroutine Package. (Input to the CORRE sub- 
routine is a linear array vice a matrix to simplify 


dimension requirements). 
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(2) The maximum number. of data cards that the program 
can handle is 1600. Capacity may be increased by chang- 
ing the DIMENSION and DATA statements. 
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A. 


APPENDIX B 


Cluster Data Program 


Purpose: The purpose of this program is to calculate 

and provide a punched card deck for input to the distance 
matrix program in Appendix C. Each input data card is 
normalized to a unit vector for the five variables NORM, 
DMHR, DML$, DOH$ and NDAY. Then each variable is stand- 


amdized by the following transformation: 


EDS _ ‚th ‚th 

where Xi; = ij— observation on the j— 
7j variable 

i. £z zn 
d» = ИШ > 
n = number of data cards 
g; = standard deviation of the ga variable 
u; = mean of the jËB variable 


Imputs: The first card ds a control card which Isa az 
lowed by a deck of aggregated data cards. The aggregated 


data cards are described in Appendix F. 


Control Card: CC 1-3 number of aggregated data cards 
Aggregated data card: the format is specified in 


Appendix F. 


Outputs: Listings of input data, normalized data and 


normalized/standardized data and a deck of punched cards 
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for input to the distance matrix program in Appendix C 


are provided. The punched output deck is in the follow- 


ing format: 
CC 
CC 
CC 
CC 
CC 
CC 
CC 
CC 
CC 
CC 


D. Special Notes: 


1-2 


4-10 
M 
Пс 
14-24 
25-35 
36-46 
47-57 
58-68 


engine/aircraft type code 
blank 

Identification code 

blank 

work code 

NORM 

DMHR E 
DML$ 

DOH$ 


NDAY 


(1) The maximum number of aggregated data cards is 150. 


This capacity may be increased by changing the appropri- 


ate DIMENSION and DATA statements. 


(2) This program may be modified to handle raw data cards 


by changing FORMAT statement number 9001. 
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APPENDIX C 


Distance Matrix Program 


Purpose: To provide a punched card deck of squared 
Euclidian distances between each pair of observation 
vectors (i.e. normalized and standardized input data 
E. The distance matrix is a required input for the 


clustering program in Appendix D. 


Inputs: A control card is required followed by a deck 
of normalized and standardized (N&S) cards. The formats 


are as follows: 


Control Card: CC 1-3 number of N&S cards 


N&S Deck: format is specified in Appendix B. 


Outputs: The outputs are listings of input data and the 
squared Euclidian distance matrix. ~Thevarsuance sia. 
is stored in a linear array that corresponds to the lower 
half of the matrix. The listing and data cards of the 
distance matrix are printed and punched in Che пва 
array format. Since a complicated procedure is required 
to locate the distance between any particular pair of 
observation vectors the format will not be specified in 
detail. Each distance matrix output card contains 10 
distances. The total number of distances calculated and 


punched is: 


n x (n-1) where n = number of N&S cards contained 
2 in the input deck. 
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D. Special Notes: 
(1) The output deck must be kept in the exact order in 
which it is punched. The position of each number in the 
linear array determines which pair of observation vectors 


were used to compute that particular distance value. 


(2) The maximum number of N&S input data cards in 200. 


This capacity may be increased by changing the DIMENSION 


statement. 
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APPENDIX D 


Clustering Program 


Lunpose: The purpose of this program із to cluster 
together individual observations and groups of observa- 
tions that are mathematically similar. Further discussion 
on clustering can be found in Ref. 8. The program begins 
by considering each input vector (observation) as a 
cluster with one observation. The two "closest" (based 

on Euclidian distance) clusters are combined into a single 
cluster and the distance from this new cluster toral 
other clusters is recomputed and the distance matrix is 
adjusted. The clustering continues until all observa- 
tions are in a single cluster. 

The program provides a detailed listing of the clus— 
ters that are combined at each step and at the users 
option the program will provide a detailed description 
of each cluster. A deck of punched cluster "descriptions" 
is produced at intervals specified by the user. A cluster 
"description" consists of the number of clusters, the 
number of observation vectors (i.e. number of elements) 
in each cluster and an ordered vector of observation 
numbers (an observation number is the number assigned to 
each observation vector based on its position when the 


distance matrix was produced). For example: 


fel 





Suppose we have 6 observation vectors clustered as 


follows 
Clusterino. No. of observation vector 
2 153 
2 2 
3 5,6939 


The cluster description would be 
number of clusters = 3 
number of observations in each cluster = 2,1,3 


vector of observation numbers = 1,3,2,5,6,4, 


B. Inputs: A control card and the distance matrix produced 


by Appendix C are required. The formats are as follows: 


Pontrol Card CC 1-3 number of observation vectors 
cc 4-6 stage at which a complete cluster 
description begins in the output 
listing. . (stage is the number of 
clusters currently existing in the 


program) 


CC 7-14 Beta value (See Ref. 8). The 
value -0.25000 used in all clus- 


tering done in this thesis. 


CC 15-17 КМ1 
СС 18-20 КМ2 See below 


CC 21-23 KM3 
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Cluster "descriptions" are punched starting at Stage КМ1 


and continuing to Stage KM2 in increments of KM3. For 


example: if KMl = 16 
KM2 = 12 
KM3 = 2 


buen cluster "descriptions" will be punched at stages 


mor, 14 mand 12. 


Outputs: A listing describing clustering as discussed in 

Bart A is provided. A punched deck of eluster deserip- 

tions aS specified by the control card is provided” me 

Punched output format is as follows: 

First card: CC 1-3 number of clusters 

Second group of cards (as many as needed): number of 
opservation vectörs In each ce lusterzunichzesehwer ee 
having a maximum of 20 numbers, each right each adjusted 
w a column field. 

Final group of cards (as many as needed): ordered vector 
of observation numbers as described in raro Ne 


format is the same as for the second group. 


Special Notes: 

(1) The maximum number of observation vectors is 100. 

(2) The user must decide at what point further clustering 
ceases to be useful. The computer algorithm automatically 


proceeds until all observations are in a single cluster. 
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APPENDIX E 


Cluster Heallocation Program 


Purpose: This program ensures that each observation 
vector is assigned to the closest group centroid calcu- 
lated from a cluster description supplied by the cluster- 
ing Program in Appendix D. The program takes the normal- 
ized and standardized data cards and based on a specified 
cluster description a group centrold is салса е а? 
each cluster. Each observation vector is then assigned 

to the closest centroid based on Euclidian distance. If 
any observation vector moves to a different cluster from 
its present assignment then the group centroids are 
recalculated and the observation vectors are again assigned 
to the closest group centroid. This reallocation process 
continues until a stable configuration is reached (i.e. no 


observation vector shifts clusters). 


Inputs: A control card is followed by Uhe normalized 
and standardized data deck then the cluster description(s). 


The formats are as follows: 


Control card: CC 1-3 number of N&S data cards 
N&S Data Deck: format is the same as specified in 
Appendix B. 


Cluster Descriptions: format is the same as specified in 


Appendix D. 


Outputs: The normalized and standardized data is printed 


along with the initial cluster description. The program 
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then prints a series of cluster centroids and cluster 
descriptions until a stable configuration is reached. 

Then a deck of reallocated cluster descriptions is punched. 
The reallocated cluster descriptions have the same format 


as specified in Appendix D. 


Special Notes: 
(1) The program is set to handle multiple cluster des- 
eriptions. Place any additional cluster descriptions 


behind the first one. -> 


(2) The maximum number of N&S data cards is 100. The 
DIMENSION statement may be changed to accomodate larger 


capacity. 


(3) The maximum number of clusters that can be specified 
in a cluster description is 50. This capacity may bë 


changed in the DIMENSION statement. 


(4) Normally about 5% or less of the observation vectors 
well Change clusters. If the percentage exceeds this 
Mt then an input error or a very unstable cubren 
configuration may have caused the problem. The results 


should be carefully checked before continuing. 
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GO BACK AND REALLOCATE ANY REMAINING CLUSTER DESCRIPTIONS 





APPENDIX F 


Aggregation Program 


Purpose: This program combines groups of raw data cards 
into aggregated data cards. A listing is provided for 

the new aggregated data. The aggregation procedure con- 
sists of calculating a numerical average for each variable 
(excluding airframe change man-hours on aircraft raw data 
cards) on the raw data card. Two additional fields are 
calculated and punched into the new cards: NDAY and NOBS 
(number of raw data cards that are combined to produce 

the aggregated data card). These aggregated data cards 
are used as input requirements for programs in Appendices 


ЕС, К апа т. 


Inputs: A series of control cards are required. Each 
control card must be the first card of thefgroupkt emac 
aggregated. The control card is punched in CC 1-3 with 
the exact number of raw data cards following it. The 
raw data cards are in the format specified in Appendix A. 
The input deck as described is called the preaggregation 
deck. It is used as an input deck for the program in 
Appendix J. Once this deck is created it should be kept 
in addition to the aggregated data cards. (Note: the 
control cards are referred to as I3 separator cards in 


the program) 
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emtputs: The output consists of a reformatted raw daca 
card with two additional fields, NDAY and NOBS. These 
new cards are referred to as aggregated data cards. 

(Each aggregated data card is an estimation of a NARF 
process.) The program provides a listing of the aggregat- 


ea data cards. The aggregated data card format foltews. 


CC 1-2 Aircraft/engine type 
CC 3 blank 
cc 4-10 identification code 
CC 11 blank 


CC 12-13 work type 
CC 14-22 NORM 

CC 23-31 DMHR 

cc 32-40 DLB$ 

cc 41-49 DML$ 

сс 50-58 рон$ 

CC 59-67 NIFR 

CC 68-76 NDAYS 

CC 77-80 NOBS 


Special Notes: 
(1) The maximum number of control cards is 150. This 
capacity may be increased by the change of the DIMENSION 


and DATA statements of the program. 


(2) The first 13 card columns of the first raw data card 
Will become an "identifier" and are transferred to the 


aggregated data card. 
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APPENDIX G 


Dominance Program for Aggregated Data 


A. Purpose: The purpose of this program is to take aggre- 
gated data cards as an input and then display in a matrix 
the dominance relationships between the data cards. 
Recall that each aggregated data card represents an 
estimated process for the NARF. The specified input/ 
output variables (in Chapter 2) are read from the input 
deck and then each process is converted to a unit output. 
The dominance matrix is printed and can be inspected for 
dominance among processes with the same output. The 
procedure to eliminate the dominance is described in 


Chapter 2. 


B. Inputs: The program requires a series of control cards 
to be inserted among the aggregated data cards. The 
control cards are punched with the number of aggregated 
data card following them (and prior to the next control 
card) in card columns 1-3. These control cards cause the 
aggregated data cards to be read in small groups thereby 
reducing the size of the dominance matrix. The control 
card should contain a number between 10 and 20 for most 
groups. For example, if we have 72 aggregated data cards 
then six control cards with 12 punched in CC 2-3 would 
suffice. However, the user must ensure that all processes 


with the same output are contained in the same control 


87 





card group. The dominance matrix is only calculated for 


those processes in the same group. 


C. Outputs: The output consists of a listing of the aggre- 
gated data group (input/output variables only) and a 
listing of the unit output processes associated with 
ch data card in the group. The dominance magix: 
printed with the above listings. Each new group is start- 
ed on a new page for ease of reference. The dominance 


matrix consists of three symbols defined as follows: 


(1) - means there is no dominance relationship between 
process A and process B (A and B are defined below) 

(2) X means that process A is dominated by (i.e. strictly 
greater than) process B 

(3) O means that process A dominates (i.e. strictly less 


than) process B. 


Each symbol is in a particular row and column of the 
dominance matrix. (the column numbers are the same as 
the row numbers since the dominance matrix is always a 
square matrix.) "A" is the row number of the symbol 


and "B" is the column number of the symbol. 


D. Special Notes: None. 
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APPENDIX H 


Listing of Aggregated Aircraft Data 


The following listing provides the values of the aggre- 
Bated data cards used in the aircraft model. The columns 
Ere as follows: 

(1) aircraft code 

(2) identification code 


(3) work code 


(4) NORM 
(5) DMHR 
(6) DLB$ 
(7) DML$ 
(8) DoH$ 
(9) NIFR 
(10) NDAY 
(11) NOBS 


The format of these cards is described in Appendix F. 
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APPENDIX I 


Listing of Aggregated Engine Data 


The following listing provides the values of the aggre- 
gated data cards used in the engine model. The columns are 


as follows: 


(1) engine code 
(2) identification code 


(3) work code 


(4) NORM 
(5) DMHR 
(6) DLB$ 
(7) DML$ 
(8) DoH$ 
(9) NIFR 
(10) NDAY 
(11) NOBS 


The format of these cards is described in Appendix F. 
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APPENDIX J 


Average Cost of Labor Program 


A. Purpose: This program computes the average cost of labor 
in dollars per man-hour for each aggregation group (i.e. 
the raw data cards being aggregated). A weighted average 
labor cost is then computed based on the number of obser- 
vations in each group. This weighted average is used in 
the Linear Economic Model to convert man-hours into a 


dollar cost in the objective function. 


Ee Eputs: The input data for this program is the preag- 
gregation deck described as an input to the Aggregation 


Program in Appendix F. 


C. Outputs: A listing of the average cost of labor for 
each aggregation group followed by the weighted average 


Est of labor for all the groups. 


D. Special Notes: None. 
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APPENDIX K 


MPS/360 Data Preparation Program 


Purpose: This program was prepared to assist the user 
in preparing input data decks for use with MPS/360. The 


required input formats are given in Chapter 4 of Ref. 7. 


Inputs: A special sequence of input data cards is re- 
mined. The input data cards must be in the order spe- 
cified. The number of data cards in each of the below 
sections is not important because the prorrateo 
cards sequentially. For example, if a particular format 
for data requires a 10 number field with 5 columns for 
each number and the input data is 15 numbers then the 
program will read 10 numbers from the first Input cara 
and 5 numbers from the second one. All remaining columns 


of the second card should be blank. 


(1) The model name is punched in CC 1-8. It must be 
alphabetic and not begin with "X". It must be left 
adjusted in the field. If right hand side (RHS) para- 
metric study vectors are required then put a "1" in CC 9. 
The RHS parametric vectors are "change columns" and are 
described in Ref. 7. They consist of a RHS column vector 
with zeros in every row except a "-1" in one of the 
resource rows. Three separate decks are created with 


one corresponding to each resource row. These decks 
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are used with MPS/360 to provide the parametric studies 


on the resource vector R. 


(2) The number of rows in the Ti matrix is punched in 
CC 1-3. The number of columns of Ti is punched in 
"om. 


(3) The penalty costs must be punched in this set of 
cards. Each card may have only 10 numbers with a 5 
column field width. The computer program automatically 
places a decimal after the first three digits in any 
field. For example, if 12282 is in the 5 column field 
then the penalty cost is $122.82. There should be one 
penalty cost for each column of Ту: 


(4) Each column of T. must be assigned a reference number 


1 
for MPS/360. Begin with "1" and continue to number 
columns until all are numbered. These numbers must be 


punched 40 numbers per card with a 2 column field width. 


ey must be left adjusted. 





(5) Each column of T, produces some particular type of 


1 
output. To associate each column with an output a set 
of cards punched with "Т, ROW" as described in Chapter 


IV must be provided. Use the same format as in (4). 


(6) The production vector (in man-hours) must be punched 
in this set of cards. The proper format is 10 numbers 
per card with a 5 column field width. There should be 


as many numbers as rows of Ti: Use only integer values. 
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(7) The resource vector, R, is specified in this card. 
Use a 10 column field width for the three values of R1, 


R2 and R3 respectively. Use only integer values. 


(8) The last input is the deck of aggregated data cards 
whose format is specified in Appendix F. They must be 


in the proper order to correspond with (3) and (5) above. 


Outputs: A deck of cards is punched for use with MPS/360. 
However, the user must provide the MPS/360 control deck 


as described in Refs. 6 and 7. 


special Notes: 
(1) Be sure not to exceed the field limits specified 


above. 


(2) The maximum number of rows and columns for the T4 
matrix is 99 in this program. See (2) above. This 

also limits other inputs to a maximum of 99 input numbers. 
Change DIMENSION and FORMAT statements if further capacity 


is desired. 


Tou 








109N*T=1 002 OG 


2 
SYOLISA ALIAILOV 3ZIWWYON 2 
ANNILNOD = 0OT 
(8I'c*ctdG'ev'ay*2zv*X2*sI*. a) LVWYOS 16 
(I)8N' (S*I=C* (I$P)X)* (I)839*(I)Əri(I)gIfI (26°59) 3119M 
(S6°9) 3LI8M (T2°03*1) dI 
JOIN‘ T=I_OOT OG 
(//* sSGONa*XL4 aSAVGNe SXTT $HOGQi XTTT 
*a$1WQs *XTT* UHWQOs *XTL* 4 WHONS ^ X85 NNOT LVOIHIINSOIaX9* Ta) LVWYOS $6 
($6*9) 3118M > 
VIVO INdNI 3LIYM 2 
(71'2*64*X6*7*632'X6'2*6J2'cv'av'zv) 1vWWu0J m 
(3109N'T-2I*(C(I)8N*(G'I- P *CI *n)X)* (1 08y*(1)88* (1)81)) (41*5) GV3% 
(O°OTJE) 1vWWuOd 21 
(£€*I-I*(1)9) (215) avay 
(0*5301) 1vWy0Y IT 
(MOYN*T=1*(I)SHE) (115) СУЗУ 
(240%) 1VWWuOd OT 
CIOON*T=14C(1)1) (OTS) Gv3% 
(240%) LVWYOS 6 
(TOONS T=I4(1)W) (64S) Gv" 
(2*G30T) 1VWY0W 8 
(102N*T=1* (1)9) (84S) УЗУ 
(£12) 1vW9u0d L 
309N*MOUN (L5) QV3u 
(2*€3) 1vWy03 S 
d (S*'S) QV3W 
(11942)  1vWyYO3 y 
VuVdW*(Z*I-I*(I)8V) (6666=0N3*9*5) QV3U 2 
VIVO INdNI OYJ% 2 
/0°T+/4NO Viva 
ӨГ 9x<1V3% 
JEW 2xu393lNI 
(Z)SV*(CE)U* (O0L )SHU^(OOT) 3* (OOT)W. NOI SN3WIG 
(001)8N*(OO )83* (OO TD) 80* (00D 81* COO T* S)X* (0T )9* (OOT C80 NOISNSWIO 
x X ok oi oc oj aic oc oic oco oic ote oc al oic o ac of ale ole aic oi ole oie ofc o oie o ade a ale к с le ad a de ade e e ol oie ok oc ooo K K SK 3k SK 3k a ad afc ok ok akc og 35 УКК с К ak 


102 








Ri ROW 
R1 ROH 


> 
° 
© 
= 
= = 
о Q = 
ac a о 
a 
u ш > 
— e 
+ >< 
- Ly 
= > N — 
O Q II - 
Ге 44 c "9 ч 
= мч 
> > = М 
со Oo u e 
O О me > 
= e”) N 
N Ne U) — 
м. x > = 
e e е >ч = Ə 
+ O O # em 1 с 
< Z Z —u+ > > 
N —-— 0 e 
- O O O en O IN N 
Га 4 XA tau кч кч om е 
N rari) Y с\) сс С\ 
mas 9 met or 
am LU ma ~ ~ a) O «~U mU 
NS eN ¿e сос Z cele 
- ОМ OON~O = te Ke XK 
CZ ome OZ ex =N I NO et 
II rt taa Tene 
i a d eo —— = == GUNnard 
eti) © oO Mx EIS e Tr 
O O — CD ez o9 со сее Фе Ф 
“O +: =c = = = — = (fe 
De ~ exea e N eja moe 
< = e — 9 +< = Nun 
~ o AE >< x *—>x<—> 
ts O cC O = О сн ч nme М D 
———— Qe wee m LL) — ees 
О О. Zae мае и Оо 6 e 
Ne De ode Or Ow JOs Os 
Nan shat AM MAIM Mm 
ee eo С || te ZH ao ce > == = 


ww “w nz u с =» Du WV) ~o 
O yb m eZ eI Fe 
A ANO RA CIs ee 
eZee EOS SR о 5 
mmx C "жес чү 
XOXO OLO coco Or OX r OO 
Zu Zu Сеи Q:xU.Ogzxuo0-zu.zu. 


= e = 
= 
o o O =O N o ~ 
n о о OO o о © 
N N “ч am сс) am 


103 


PUNCH INPUT DECK FOR MPS/360 


= 
~ со 
- O 
+ - 
<= e 
N > 
Oe N — ~ == 
Orm e N rt —ч 
н = < = ~ 
ft = e e = 
za e«t - - - eo 
MAN = L m O 
Dr e > - - e 
zu чі — = - ~ m~ 
NGO e X >< - =ч 
te NO MN e UV) € 
Su. oze мч e = J 
Y = eO e е sO 
LUNG — —J =) DUE 
OZOZOZZ e Nee OZ 
"NINO eQ = eC e WO =O 
СУМ сб суо e fx Ou» NOOO 
mo oe ee || om II e— oo u e. 


mu m um we mu gem ul mus m = 
w w “wx — D ~ Dw ~ 
e -e eO r-ZO eZ FO 
W <T U <1 U < O U A e OU IA LIA Ош 
be SR SPR EEN SE EF ZO 
mate Ome ке € 
e а: о ceOcdccocla c с 
ZU. ZLZLOZLOOZUL-VOUELOS 


O O O 0 AOD O 
= су Ln оо оо N 
e m m re AN Ln 





“=ч 
t 
° = 
N + 
e ° 
ц. < 
e et 
x u. 
со - 
e >< 
- со 
LE e 
Га 4 - 
- e =~ 
e (Y — = - 
>< - - - O 
e e © O e 
e >< e e e 
Es dm ~ ~ o o L 
+ ° - ~ + = - - 
eN + o e > e e 
Ne © N N x x x 
=ч *Ш (y = e m + m 
UN entr ц. и. = = oa 
e >< MLL - - - - - 
XXI e > x N ~ eo 
= = =—>< m со < =з =з 
eL WN >< CO e e e e e 
- =< = e N =] - - - 
co ee N еч & es + - - - 
C LI! eco -— - о - e е 
e. Ne „ш X x x x 
un Te = O Oo O vo O 
> = >x >x >< e * ~ e j ~ =e + e e 
m om € o> MX met oO rm 
or Hor Фф ~ М wf O w | — pt 
NN С\ me C = «5 -j ~ Ом е Mea = 
{ue <í < <T w. е е es E ro oe - 
end onm е ж ч. N) m EIN IN Im < 
- = 9 On т TO e з= O I T = — 
NN-N mow Yung ~ eo Caw BOX agx aq O 
e Ne Me 2 OwZwurme M~e OmZrie M mie se OO OA 
eQ eO e NT WO = eOe We ete ene fe WZ Oe 


XOXO Њиыгом Ох et HOOK AY OK ORM ОШ —— 


$ eF eF se Zef I eF аи оф иж фо еф ce “š 
—— P DD Sep NS шна Мешок ә S — UL! sZ Cy >< cy Qv 


e e eZ Oana EZO r2000 20 • SL EZ o 

a o A AS O ш 
ZELEEEEBZO-LREEOLEZEeOQrLEZEQOo-6rEerLEEBZELELEZrLEZaA 
Lm mX Zm Se — = mY MY SZ OOO 
Odc Oo OOdcOOcIcOOOdcOOu. OQOdcoootuowoOoOc«oooccoocdcor-z 
WSL SUL (OE ULO QC 7E LL OO O 7S UL (0 - OO zx UL CO Of XUL O ZU OO 7X uL Co o Uu. LU 


an 
c4 N mo O HO HO =O 4 OAODOONO OO 
Ө О со © OO oo NN M TAATA po 
© о 00 N 00 oo OO о СС ОООО Or 


104 





APPENDIX L 


Cluster Centroid Dominance Program 


A. Purpose: This program takes an input deck of aggregated 
data cards (processes) and a set (more than one if de- 
sired) of cluster descriptions and produces a dominance 


matrix showing the relationships between cluster centroids 


(or "aggregated" processes). 


B. Inputs: A control card with the number of aggregated 
data cards punched in CC 1-3 is first. Aggregated data 
cards then follow this card. Next the user may put 
several cluster "descriptions" if desired. (At least 
one is required.) The cluster "description" format is 


specified in Appendix E. 


C. Outputs: A cluster description, centroid data listing 
and dominance matrix are provided for each cluster des- 
eraption. A list of efficient groups is ее M 
efficient group is a centroid that is not dominated by 


any other centroid. 


D. Special Notes: 
(1) The dominance matrix shows an "O" where each centroid 
is either dominates or is dominated. The number of "O's" 
in each row and column are listed at the sides of the 


matrix. A "-" means no dominance relationship exists. 
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. (2) The maximum number of centroids in any cluster des- 


cription 1s 50% 


(3) The maximum number of aggregated data cards is 100. 
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